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1 Introduction 

A full understanding of the infrared structure of QCD matrix elements is an unavoid¬ 
able step towards making precise predictions of Standard Model backgrounds at hadron 
colliders. In order to make finite predictions for cross-sections we must cancel infrared 
singularities between unresolved real radiation corrections and singularities in the virtual 
(loop) corrections. The study of infrared properties of perturbative gauge theories have 
a broader scope beyond this application since the universal behaviour provides a strong 
constraint on the structure of scattering amplitudes. 

The soft and collinear infrared limits at next-to-leading order (NLO) have been under¬ 
stood long ago and general algorithms (e.g. Catani-Seymour [1] or Frixone-Kunszt-Signer 
[2]) for the computation of infrared finite cross-sections form the core of the current gen¬ 
eration of precision tools used to make theoretical predictions for the LHC experiments. 

In the last ten years or so a lot of effort has been put into generalising these techniques 
to next-to-next-leading order (NNLO) and a variety of different techniques now exist with 
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the ability to make finite predictions for important LHC observables (e.g. references [3-8]). 
All of these methods rely on knowledge of the underlying factorisation properties of QCD 
amplitudes in the double unresolved limits at tree-level [9] and single unresolved limits at 
one-loop [10-14]. 

The first step at next-to-next-to-next-to-leading order (N^LO) has been taken recently 
through the complete calculation of fully inclusive Higgs production at hadron colliders up 
to 0(ag) in the large top quark mass limit [15]. This calculation has been performed in a 
number of different stages building expansions around the soft limit [16-18] and using the 
reverse unitarity method to obtain each component of the triple-virtual [19-21], squared 
real-virtual [22, 23], double-virtual-real [24-27], double-real-virtual [28, 29] and triple-real 
radiation [30] as an expansion in the dimensional regularisation parameter. The poles of 
these separate contributions cancel analytically when summed together and combined with 
the counter-terms for UV poles [31-34] and initial state infrared singularities [35-39]. 

Further steps are required to extend these techniques to fully differential observables in 
an analogous way to the NLO and NNLO cases. Many of the infrared regions that must be 
accounted for in such a procedure are now fully understood. The missing ingredients that 
remain are the one-loop triple collinear splitting functions involving gluons. Though the 
factorisation of the squared matrix elements are sufficient for the construction of infrared 
finite cross sections, factorisation at the amplitude level [40] can yield much more com¬ 
pact expressions leading to a more efficient construction of the factorised squared matrix 
element, especially when considering spin-correlations. 

Figure 1 shows the real and virtual contributions to a cross section up to N^LO and the 
primary singular limits which are either multiple soft, or multiple collinear, 

The factorisation amplitudes have been computed in all cases [9, 24, 26, 41-46] except for 
the triple-collinear and double-soft limits of the double-real-virtual. The triple collinear 
limit at one-loop has been considered at the squared amplitude level for q —)■ qQQ [47] and 
for the mixed QCD-I-QED cases of g —)• 977 , q —)• qg'^, g —)• qqj, 7 —)• qq'y and 7 —)• qqg 
[48, 49]. 

In this article we compute the one-loop gluon initiated splitting functions in QCD, 
g —)• ggg and g —)• qqg. Splitting amplitudes valid in four dimensions are extracted from 
the known analytic amplitudes for 0 —)• Hgggg [50-54] and 0 —)• Hqqgg [50, 55, 56] in the 
large top-mass limit. 

The structure of the article is as follows. We first introduce the notation for the am¬ 
plitudes and the squared amplitudes together with their respective colour decompositions 
and collinear limits. In section 3 we describe a parametrisation of the multi-collinear limit 
using spinor-helicity variables which we will use to compute the splitting amplitudes. We 
then present the g —)• ggg and g —>• qqg splitting amplitudes and describe the symmetries 
and super-symmetric decompositions used to obtain a compact representation. We then 
check the universality of the new splitting amplitudes by taking a numerical limit of the 
gg —7- gggg and gg — ?• qqgg in NJet before reaching our conclusions. 
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Figure 1: The contributions to perturbative cross sections up to N^LO. This consists of 
virtual (V) corrections up to three loops and real radiation (R) corrections with up to 3 
additional unresolved legs. In the real radiation contributions the primary infrared limits 
of soft (S) and collinear (C) should be removed from the matrix elements and re-combined 
with the virtual corrections to obtain an infrared finite result. 

2 Notation 

A general QCD amplitude can be decomposed into a basis of SU{Nc) colour factors and 
ordered partial amplitudes which depend only on the momenta and helicities of the external 
legs. For an n-point L-loop amplitude this can be represented as, 

{??■}) = E c=({“.))^S({p1‘}) (2.1) 

c 

where Oi,Ai and pi are colour indices (adjoint or fundamental), helicity and momenta of 
the M leg. Unless explicitly indicated otherwise, we understand that the index i runs from 
1 to n, e.g. 

= ( 2 . 2 ) 

For cross-section computations we are required to square these amplitudes and sum over 
the colour indices. This sum can be represented as, 

= E (A'->({»<).{p 1‘)))'xl^'>({a^),{p^}) 
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Figure 2: Factorisation of tree and one-loop amplitudes in the multi-collinear limit. 


where the matrix ^ is a function of defined by 

^ {C,({<n}))'c,.({a.)). (2.4) 

di 

while Ad) is a vector of partial amplitudes 

= (2.5) 

r^i 

Partial amplitudes may in turn be written in terms of primitive amplitudes Ap ' which 
further decompose colour and flavour structure due to the internal loops, 


Ad} = ^ R,,p,x{Nc, Nf)Adf\ ( 2 . 6 ) 

p,X 

where X runs over the independent primitive topologies at L loops and p runs over per¬ 
mutations of the n external legs. Eq. (2.3) can thus be equivalently written as 



where is a vector of primitive amplitudes A^n,^^ and the matrix ^ can be related 
to Cn^’^ ^ defined in Eq. (2.4) by the change of basis in Eq. (2.6). 

In the limit where m of the external legs become simultaneously collinear, the am¬ 
plitudes factorise into a product of lower multiplicity amplitudes and splitting amplitudes 
which contain all the infrared divergences: 


a<P({pI-}) ‘"4 


L 

EE' 

/c=0 Ap 


SpS-'-'t-P-"'-; {pi'<)”i)44.+l(^’"'’. (2.8) 
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where and can either be primitive or partial n-point amplitudes and splitting 

amplitudes respectively, while and P = pi + ■ ■ ■ + ■ A schematic representation of this 

factorisation is shown in Figure 2. The sum of internal helicity states Ap leads to spin 
correlations in the factorized squared amplitude 


L L' 


k=0k'=QXp,X'p 

where we can define 

{pi'}))* ■ ^<'■•'■'1 ■ 


T^i.ppp-.ipp) = (s1>i'''(P"';{pf 

in terms of partial amplitudes or equivalently 


-i{L,L') ) p dAe>/. (.Pi 


■'Sp,n 


Sp; 


= (sp^V"";{ 4 }))^ -cgf • ^Vn\p^^'\{p-A) 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


\L L^l 

in terms of primitive amplitudes. In the colour matrix Cgp’^ we absorbed a prefactor 
which takes into account colour conservation along the factorized parton, such that 


niPP] 


1 


1 V 2 - 1 
1 p[L,L' 


r\PP\ 

On 


for gluon-initiated Sp 
for quark-initiated Sp 


(2.14) 


and similar for Cgp’^ ^. 

For brevity, the results presented in this paper will often omit the subscript indicating 
the number of partons involved in an amplitude, since this can be deduced by its arguments, 

i.e. 


Sp(-P^";4,... ,4") = Spj-p"- fe"“}) 

A(p^,...,p^) = An({p,^*}). 


(2.15) 


3 A spinor parametrisation of the multi-collinear limit 

We define the multiple collinear limit using a parametrisation of the full kinematics in term 
of a parameter 6, such that the collinear limit in Eq. (2.8) is identified as the leading term 
as (5 —>• 0, i.e. 

lim = lim4^n{4(<^)}) 

o^U 

= EEsp!»-*>(-^’-"mpl)£i)4'‘-™+i(p"'.{pl)?™+i)+o(^). (3,1) 

k=0 Ap ^ 
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The parametrisation is defined by, 


P>t{5) = ZiP^ + 6k!^ . - 6^ 


^T,i 


2{P ■ r])zi 




Pii^) = Kl'id, {p}m+l,n,v) 


i = 1,... ,m (3.2) 

i = m + 1,... ,n (3.3) 


where Zi = {pi ■ p)/{P ■ p) are the momentum fractions of the unresolved partons, p is an 
arbitrary light-like momentum and P is the massless projection of P = 

P^ .. 


pfJ. _ pM _ 


-p^ 


2P-P 

The vectors are orthogonal to P, P and p 

kT,i • P = kT,i • P = kT,i - p = 0. 

Momentum conservation implies that: 

m 

^zi = i 
i=l 
m 

E 4 ,. = o'‘ 

i=l 

m 


'^T,i 


= P\ 


(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 


i=l 


The function is a generic map that keeps the factorized momenta m -|- 1,..., n on-shell 
as well as absorbing the recoil P^/(2P ■p)p^, and it satisfies pf as (5 —)• 0. The exact 

form is not important for our purpose of explicitly taking the limit and various mappings 
have been considered in the literature (for example in the Catani-Seymour subtraction [1] 
or Kosower’s antenna [57]). When implementing the collinear phase-space numerically we 
employed the Catani-Seymour map as described in Appendix A. 

Since we are working at the amplitude level, we would like to have a parametrisation 
of the limit valid for the spinors of pi as well. This can be achieved using an appropriate 
choice of the transverse vectors kT^i, 


26k!^ - = {zi)[u}i]{P\j>^\p] + [zi]{uJi){p\-f>^\P]. 

In the above we use the notation 


(zi) = 


(ip) 

{PpY 


Zi = 


[ip] 

[PpY 


{^i) = 


jiP) 

{ppy 


m = 


M 

[pP] 


(3.9) 


(3.10) 


where the spinor variables {zi) and \zi\ differ by a phase from the usual parametrisation 
which uses y^. It is worth to notice that both {oji) and [wj] are 0(S) in the collinear limit. 
The spinors parametrisation then reads. 


|i) = (zi}lP} + (ujYIp) 


|i] = [zi]\P] + [uji]\p\. 


(3.11) 
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We find that this is a convenient way to take the limit at the amplitude level since the 
spinor variables {zi) obey Schouten identities: 

(z,)(jA:)=0, (3.12) 

ijk cyclic 

as well as momentum conservation, 

'L l>t - f“‘- ■ ( 3 - 13 ) 

i 

For the triple collinear splitting amplitudes this means we have the kinematics of a five- 
point function event though the colour space is that of a four-point function. 


3.1 Example: the tree-level MHV multi-collinear splitting amplitude 

The result for the multi-collinear limit of the maximal-helicity-violating (MHV) amplitude 
has been known for a long time. More recently the general helicity cases were also examined 
through use of the MHV rules [45, 46]. This case is incredibly straightforward and serves 
as a useful example of the general treatment introduced in the previous section. 

We start with the Parke-Taylor MHV amplitude with particles 1 and r > m having 
negative helicities and all others positive helicity. 


(l-,2+,3+,...,r-,...,n+) = 


(Ir) 


niLi(n + l)’ 


(3.14) 


where the product in the denominator is considered modulo n. The limit is simply taken 
by applying eq. (3.11) 


HW(l-,2+,3^ 


,n' = 


({zi){Pi^rar) + (wi)(r?r)) 

]Xj=i{jj + 1) {{zi){nPi^m) + (wi)(nr/)^ {{zm){Pi,mm + 1) + {ujm){vm + 1)) 0”=^ + 1) 

_ (*^ 1 ) _ (.Pl,m^) _ 

{Zm) nT=lijj + 1) {nPl,m){Pl,mm + 1) + 1) 

= Sp(0)(-P+^; 1-, 2+,..., m+) (m + 1)+,..., r",..., n+) + 


(3.15) 


where we have used eq. (3.9) to perform the power counting. For i^j G [l,m] this can be 
seen explicitly. 


(wi) 


^ (u) 


^ {Pl,m\kT\ri\ 

2(Pi,m • 


0{5) 


i{Zi){u}j) - {Zj){uJi)) {Pl,mV) 


0 { 6 ). 


(3.16) 


One can clearly arrive at this final result without being so explicit about the parametrisa- 
tion, yet it is convenient to have one in a generic implementation. 
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3.2 One-loop basis functions for pp ^ H 2j in the triple collinear limit 


The analytic H + 4 parton amplitudes have been computed using unitarity cuts and ex¬ 
pressed in terms of the universal infrared poles plus finite logarithmic and di-logarithmic 
functions as well as rational terms. Taking the triple collinear limit of the infrared poles, 
rational terms and logarithms as above presents no difficulties. Dealing with the di- 
logarithmic parts requires some minor effort to ensure the arguments are in the appropriate 
region so the limit will converge. Polylogarithmic identities are well known and understood 
in huge detail (see Ref. [58] for a recent review) - way beyond the simple structures ap¬ 
pearing here. Nevertheless we collect some potentially useful identities here to aid the 
reader, 


TT 


Li 2 (1 - x) -h LD (x) log(x) log(l - x) - — = 0 

D 

Li 2 (x) -h Li 2 Q j -b ^ log(-x)^ ^ = 0 


X G [0,1] 
X < 0 


Li2 


xy 


(1 -x)(l -y) 

X 


- Li2 - 


1 — X 


- Li2 - 


1 -y 


+ 


— Li 2 


1 -y 


— Li2 


1 — X 


- log^ 


1 — X 


= 0 


x,y G [0,1] 


(3.17) 

(3.18) 


(3.19) 


One function requiring a bit more thought is the three mass triangle which has square roots 
appearing in the arguments of the di-logarithms [59-63]: 


ir ^ (1 - - Li2 (l - i) 

- ^ log^ (Zk) - log log (Zk)^ (3.20) 

One other minor issue with the results available in the literature is that the NMHV expres¬ 
sions have been presented using Forde’s method for triple cuts [64]. This method gives the 
coefficients as the sum over solutions to the on-shell equations. To aid our computation we 
performed this sum explicitly to write the coefficients in terms of the usual spinor products 
of the external momenta. 


4 y —)■ ggg splitting amplitudes 

4.1 Colour structure and primitive amplitude decomposition 

In the section we will suppress all helicity superscripts and the function arguments are 
taken to represent both momenta and helicity. The tree-level colour decomposition can be 
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written as 


Sp^^Xiap, ai, 02, as}, -P; 1, 2, 3) 

= X] (4.1) 

o-eSa 

= ri“'^(2)6/6M3)“PSp(°)(-P;l,a(2),a(3)) (4.2) 

(tSS2 

where tr(oi,..., an) = in terms of the fundamental generators of SU{Nc) 

and 7“^'= = ^^/2/“'>'= in terms of the adjoint structure constants. The relation between the 
two representations can be shown to hold using the Kleiss-Kuijf relations [65] for the 
splitting amplitudes, 

SpW(-P; 3, 2,1) = SpW(-P; 1, 2, 3) (4.3) 

Sp(0) (-P; 1,3,2) = -SpW (-P; 1,2,3)- Sp(o) (-P; 1,3,2) (4.4) 

The one-loop colour decomposition is^, 

Sp^^\{ap, ai, 02, Os}, -P; 1, 2, 3) 

= X] tr(ap, 0^(1), 0^(2), a^(3))Sp^^^(-P;(T(l),o(2),o(3)) 

o-eSs 

+ tr(ap,aa(i)tr(a,7(2),a<7(3))Sp3^^(-P;'7(l),cr(2),cr(3)) (4.5) 

o-eSs/Za 

where the partial amplitudes are composed of primitive amplitudes as follows: 
Sp«(-P;l,2,3) 

= IVcSpt^](-P; 1,2, 3) - iV/Sp[^l (-P; 1,2,3), (4.6) 

Sp(')(-P;l,2,3) 

= 2 (Sp[^l(-P; 1, 2, 3) + Sp[3l(-P; 1, 3, 2) + Sp[3](-P; 3,1, 2)) . (4.7) 

The primitive amplitudes for the gluon and fermion loops obey line-reversal symmetry, 

SpW (-P; 1,2,3) = Spl^l (-P; 3, 2,1) (4.8) 

and so in all we have three independent gluon loop primitive amplitudes, three fermion 
loop primitive amplitudes and two tree-level primitive amplitudes. The colour summed 

^We write the one-loop decomposition in the standard trace basis rather than the slightly more compact 
‘F-basis’ representation of Del Duca-Maltoni-Dixon [66]. Since we express the colonr summed squared 
matrix element in terms of the minimal basis of primitive amplitudes the final expressions are eqnivalent 
to the DDM forms. 
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Born and virtual corrections can then be written according to (2.13) using: 


-^[0] _ /SpM(-P;l,2,3)\ 
l^SpM(-P;l,3,2)J 





/iV,Sp[3l(-P;l,2,3)\ 

iV,Sp[3l(-P;2,l,3) 

iV,Sp[3l(-P;2,3,l) 

iV/SpW(-P;l,2,3) 

iV^SpW(-P;2,l,3) 

ViV/SpW(-P;2,3,l)/ 


^[ 0 , 1 ] 

•-Sp 


2 iVc 


/2 -2 0 -2 2 0 
VO -2 2 0 2 -2 


(4.9) 

(4.10) 


(4.11) 


(4.12) 


We also choose to present the results using the super-symmetric decomposition: 
Sp[^l(-^’;1,2,3) 

= Sp[^=^] (-P; 1,2,3)+ 4Sp[-^=0 (-P; 1, 2, 3) + (1 - e5ij)Sp[^=°l (-P; 1,2,3) (4.13) 

Sp[^l(-^;1,2,3) 

= Spl-^=0(-P;l,2,3) + Sp[-^=°l(-i^;l,2,3) (4.14) 


since this yields particularly compact expressions. We also include the scheme dependence 
for both the FDH = 0) and CDR {Sr = 1) schemes. 


4.2 Results 


We define the following phase-free quantities. 


C^ij — C^ijk 


{ij){zk) 

{jk){zi)' 


Pij — Pijk 


[jk][uJi] ’ 


[+1 


(4.15) 


Since there can be no repeated index in either aijk and fiijk each can be uniquely specified 
by the two first labels. 
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The integral functions are defined using the following basis, 


^ ^ (^ 3 ) + y ) - log 
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2 

+ log 


vr 
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S 123 


log 


S23 

S 123 


1-^3 

Zl 


log 


2:2 


2 l 


+ Li2 I-) + Li2 I-) + Li2 I — 


■§12 

■5123 

22 


+ log 


1-21 


23 


log 


■§23 

•5123 


23 


23 


— Li 2 1 — 


Sl2 


-V 


(1 - Z 3 ) 5123/ 


Li 2 1 — 


1-23 

S23 


+ Li 2 1 — 


2 l 


1-21 


(1 - Zl) 5123 ) 


Fr'"' = -los(l-z,) (log( 


^1^3 A 'Sl2 \ 


- - ( log (zs) log 


1-23 

^ + log ( 21 ) log 


+ log y—j ) + log (^ 1 ^ 3 ) log 


S 12 

S 123 


77NMHV ttNMHVI 
^2 = 


lo3 


R 


3™V = l(log(z3)log 


•5123 


■§12 

•5123 


■§23 

Sl23 


TT 

~3 


+ log( 2 ;i)log 


- log 


■512 

Sl23 


)‘° KS ) 


^bo” = ---logM^)-2(Li2(l- 


TT 


2 / 'S 12 


5123 

S 12 


523 \ ^ 

5123/ 3 


+ Li 2 ( 1 — 


5123 \ 
523 / 


Lo ( 51 , 52 ) = log 

Li ( 51 , 52 ) = —^— log ^ 
51-52 \S2j 


log ( - ) - ^ ^ 


1 1 

9 1 - 1 - 

Z 5l — 52 V 5l 52 


L 3 ( 51 , 52 ) = 


log ( 


5l 


(51-52)^ V' 52 / 2(51-52)2 Vsi '52 


1 1 

— + — 


(4.16) 

) 

(4.17) 

(4.18) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 


We express the infrared poles and associated logarithms as described by Catani’s for¬ 
mula [47], 

•4 = -K(^)‘ + (^)‘ + ^ 3 -' - 2 )) 0 - 2 ^) 

All results in this section are presented unrenormalized. 
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The tree-level splitting amplitudes are 
Sp(o) (-P+;l+,2+,3+) = 0 
Sp(0) (-P+;l-,2-,3-) = 

Sp(0) (-P+;l+,2+,3") = 

Sp(0) (-P+;l+,2-,3+) = 


[zi] [zs] [12] [23] 
(zi)(12)(23) 


(zi)(z3)(12)(23) 

Sp<») (-P- 1+. 2-.3-) = - IT1! ( fe + 

[23]^ lysm (1 - ^ 3 ) 2 ^ 3512 / 


(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 


Sp(0) (-P+;l-,2+,3-) = - 


[2Pf 


[13]^SipS3P 


^21^2343 732^2^:3'S1P («13«3i) 7l2^1^2S3P (ai3a3l) ^ 


■Sl2'Sl23'S23 


723 (1 - 5 ^ 3 ) 


721 (1 - 2 : 1 ) 


(4.31) 


All other helicity configurations are given via parity or the line-reversal symmetry of eq. 
(4.8). The one-loop splitting primitive amplitudes are, 


Sp[-^=4] (_p+. l+, 2+, 3+) = 0 

Sp[-^=4] (-P+;l-,2-,3“) =Sp(°) (-P+;l“,2-,3“) (Fg + F^^^) 
Sp[-^= 4 ] (_p+; 1+^ 2 +, 3-) = Sp(°) (-P+; 1 +, 2 +, 3-) {Vg + F^hv^ 
Sp[-^=4] (-F+;l+,2-,3+) = Sp(°) (-F+;l+,2-,3+) (Fg + F^^^) 


(4.32) 

(4.33) 

(4.34) 

(4.35) 


Sp[A^=4] (-F+;l+,2-,3-) =Sp(°) (-F+;1+,2-,3-)(F3)- 


[23] 2 




^ 723 « 12 / 1 i 2 ' S 23 ^ ■§12 (^^1 - 1 ) ^ ^ 1 


2^3512 


/ ^12'! 


23 


Z 2 SIP 
1 3 


^^3 — 1 / 'S 123 VS12S1PS3P 


- /532 


7 NMHV 




723ai2^12'S23 (21 - 1) ^ 77 NMHV 
^2 

Z 2 Z 3 S 12 SIP 

1 f /ll2sf 


’23 


S 123 VS 12 SIPS 3 P 


+ 1^32 F 3 


pNMHV 


) 


(4.36) 


Sp[A^=4] (_p+; 1-, 2+, 3 -) = Sp(°) (-F+; 1-, 2+, 3 -) (Vg) - 


[13]2 


-h 

-h 


^NMHV 


/ 732/^31'Sl37f2 ^22:3(«13«3l)^7|2 7f2^2/?|l7|2 

V 723^3'S12'S3P 723'S3p(^ 3 “ 1) 723^3'S3p(^3 “ 1) 

(ai3a3l) ^7 i2 , I 12 P 13 P 31 I 32 , 712^1^ (ai3a3l) ^32^ pNMHV 

Po 


Sip 


-h 




723^;1'S23(5;1 - 1) Sip(^l “ 1) 

A \ \ 

iNMHV I 


S 123 \S12S23'S1PS3P / 


) 


(4.37) 
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Sp[A^=i] (_p+; 1+, 2+, 3+) = 0 
Sp[A^=i] (-P+;l-,2-,3") = 0 

Sp[A^=i] (-P+; 1+, 2+, 3-) = Sp(°) (-P+; 1+, 2+, 3-) a32Si2^i (s23, S123) 
Sp[A^=i] (-P+;l+,2-,3+) = (-P+;l+,2-,3+) 


■^box 


— -£1 ('S 12 , 5123 ) +-^1 (^23) ■S 123 ) 513 


[li^]='s23 


Sp[A^=l] (-P+; 1+, 2-, 3 -) = - (Si2, S 123 ) 

Sp[A^=l] (_P+;1-,2+,3-) =-: 


[ 23 ] 2 si 23 

[2P]2 


7 ?lm 

'^box 


[13]25123 V 
— ( Li (si2, S123) + Li (s 23 , ■S123)') 'S13 


(4.38) 

(4.39) 

(4.40) 


(4.41) 

(4.42) 


(4.43) 


SplA^-ol (-pM + .2+,3+) = (J- - (4.44) 

3(12)(23) Vsi23 P31S12S23 ) 

Sp[-^=oi (-P+;l-,2-,3-) = ^Sp(°) (-P+;l-,2-,3-) 

Zx (1 - 2:2^) 232:2 , , 2123 ( 73221 71223 S13 

2 i 22 + T--r—Tj- 7 - 7 --r + 2322 + 2123 -( -h -h 

(1 - 21) (1 - 22) (1 - 23) 


712732 V 1 - 23 1-21 •5123 


Sp[Ar=0] ('_p+. 2+, 3") = Sp(°) (-P+; 1+, 2+, 3") 


37 : 


32 


+ 


S12 V 23 23-1 

1 /S23 1 


732 2 712 -j _ ^^27321/2 (S23, 5123) ■S23 + 2712/323-^3 ('S23, 'S123) Sl3'S23 


Sl23 V'S12 2 


- 7 (712 + 2) 732 


Sp[-^=01 (-P+;l+,2-,3+) = Sp(°) (-P+;l+,2-,3+) [ 

30^31 \ 


2-//3 ('S12, 'S123) Sl3 


0:32 



5/2 

1 

+ 


H-- 


2 V^123 

•523 

+ 

•5123 

oplm 

■J-^box 


— 20)32-^3 (s 23, SI23) -513 — 31/2 (S23, 'SI23) (2si2 + 3si3) 


1 f S 12 , S 23 


2si23 V'S23 ■S12 


Sl2'S23 Sl3 


(4.45) 


(4.46) 


(4.47) 
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5 g ^ qqg splitting amplitudes 

5.1 Colour structure and primitive amplitude decomposition 

The colour structure of the tree-level splitting amplitudes is 

5p(°^({ap,n,i2,a3}, --P, Ig, 2^,3) = 

r(ap,a3):^Sp(0)(-P;l,-,2„3) + T (as, ap)| Sp(o)(-P; 2^, 1,-, 3) (5.1) 

where T(ai,..., an)l = Note that charge conjugation symmetry allows 

us to write Sp^'^^(—P; 2^, Ig, 3) = Sp^°^(—P; 2g, Ig, 3). At one-loop we have three colour 
structures, 

Sp^^\{apJi,i2,a3},-P;lq,2q,3) = Nc Sp^^|(-P; Ig, 2g, 3) 

+ r(a3,ap))i Spg(-P;2g,lg-,3)' 

+ Spg(-P;lg-,2g,3). (5.2) 
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3 


Sp[«>(-P; l,-, 2 „ 3 ) = 




3 , 



Figure 3: Sample diagrams corresponding to primitive amplitudes for Sp^^^(—P; Ig, 2g, 3g) 


The partial amplitudes Sp 4 .]^ and SP 4.3 are given in terms of the primitive amplitudes 

SP4.4(-P; Ig, 2g, 3) = Sp[^5 (-P; Ig-, 2g, 3) - ^Sp[^l(-P; Ig-, 2g, 3) 

+ ^SpW(-P;lg-,2g,3) (5.3) 

SP4,3(-P; Ig-, 2g, 3) = Sp[^l (-P; Ig-, 2g, 3) + Spt^l (-P; 2g-, Ig, 3) + Sp[^ + ^1 (-P; Ig, 3, 2g) 

+ Sp[«l (-P; Ig, 2g, 3) + Sp[^] (-P; 2g-, Ig, 3) (5.4) 

where the indices [L] and [R] label the primitive amplitudes corresponding to fermion lines 
turning left or right upon entering the loop and [/] denotes the primitive amplitudes with 
fermion-loop contribute. The label [L + P] in the sub-leading colour amplitude corre¬ 
sponds to the sum of the left and right primitive amplitudes for the non-adjacent fermion 
conhguration. Some representative diagrams of the primitive amplitudes are depicted in 
fig. 3. 

The colour summed Born and virtual corrections can be written as in Eq. (2.11), where 
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gp[scalar](_p. 2„ 3) 


+ . . . 



Figure 4: Sample diagrams corresponding to scalar contribution for Sp^^^(—P; Ig, 2q, 3^) 
in AA = 4 super-symmetric Yang-Mills theory. 


the vectors 



and the colour matrices C 


iL,L') 

Sp 


are now given by 



^(0,0) 

^Sp 



C 


( 0 , 1 ) 

Sp 


( Sp(0)(-P;l,2,3) \ 

l^_Sp(0)(-P;2,l,3)J 

-1 \ 

/ Y,Sp4.i(-P;l,2,3) \ 
—A^c Sp4.4(—P; 2,1, 3) I 
\ Sp4.3(-P;l,2,3) J 

-1 Nc\ 
^c\ -I N^-IN,) 


(5.5) 

(5.6) 

(5.7) 

(5.8) 

(5.9) 


The quark primitive splitting amplitudes also have a useful super-symmetric decom¬ 
position [11]. In this case we can write the complicated “left-moving” amplitudes in terms 
of simpler ones built using the A7 = 4 super-multiplet, 


Sp[^] (-P; Ig, 2g, 3) = Sp[^=^l (-P; Ig-, 2g, 3) - Spt^l (-P; Ig-, 2g, 3) 

- Sp[7](-P; l^-,2g,3) - Sp[«^‘^i'^'-l(-P;lg-,2g,3) (5.10) 

SpT+^](-P;l^-,2,3g) = Sp[^=4l(-P;lg-,2,3g) -Sp[^^"i'^'-](-P;lg-,2,3g) (5.11) 


where the function indicates the contribution from the complete scalar sector of 

A7 = 4. This includes a scalar-fermion-fermion vertex as well as the scalar-gluon-gluon 
vertex which contributes to the function in the pure gluonic case. Representative 

diagrams contributing to are shown in fig. 4. Using these relations we find a 

compact form for the colour dressed splitting amplitudes in terms of {~P] 1?) 2g, 3), 

Sp[^l(-P;lg-,2g,3), SpW(-P;lg-,2g,3), Sp[-‘^i-l(-P;lg-,2g,3), Sp[^="](-P;lg-,2g,3) and 
Sp[-ai-](_p. i_^2,3g). 
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5.2 Results 


As before all results in this section are presented unrenormalized. The non-vanishing 
independent tree-level splitting amplitudes g —)• qqg are 


Sp(0) (-P+;i+,2-,3-) =- 


{z2){z3) 


721 


/3i2[23](12) 1^(1 —^ 3 ) 2:3 732'Si2372iy^ 


•512 \ 


(5.12) 


Sp( 0 )(_p+.i-^ 2 + 3 -) =- 


(^ 3 ) (/33i) 
[13] (2P) 


t 


S23 (a?3/5i3)^ zf ^ S13 (731)^ ^ zisi2 


Z3Sl3{P3lV 


Z 3 S 123 (1 - Z 3 ) Si3 


(5.13) 


Sp(0) (-P+;1+,2-,3+) = 


{Z2r 


{z3){12){23) 


SriW ('-P+-1" 2+ 3+) = 

P [ ^ ) (^ 3 )( 12 )(23) 

and the others are obtained by conjugation using the relation 


Sp(-P+;lJ^2^^3'^3^ =Sp(-P-;l 


— 9—^*2 Q—hs 

q 1 1 




(5.14) 

(5.15) 

(5.16) 


The sub-leading colour tree-level splitting amplitudes g —)• qgq are not independent because 
they can be expressed in terms of (5.15) using the KK relation (4.4) re-written with the 
quark labels, 


Sp(0) (-P+;lJl,2'*^3^3) = 

-Sp® (-P+;lJ\ 3^3^ 2^*2^ +Sp(°) (-P+;3J^1^^2^2^ (5.17) 

A sample of two representative tree-level splitting amplitudes g —)• qgq is 

Sp<«)(-P+;l+2-,3-) = (5.18) 

Sp<«)(-P+;lP,2+,3l) = -^AT_ (5,19) 
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The non-zero independent one-loop splitting amplitudes g —)■ qqg are 


Sp[A^=4] (_p+; 1+ 2-, 3-) = Sp(0) (-P+; 1±, 2-, 3-) 

_ {Z2){z‘i) ( (1 - Zif Si 2 pNMHV _ f ihslp ( 712 )^ _ T 23 \ pNMHV 

[23](12) y 7215;3'S1P ^ y (1 — 2 : 3 ) 723S12'S3P ('i- — Z3) Z3 J ^ 

_ ^ 3512^13 _ / Zil 3 i 2 sl 3 _ 7 |lglp( 73 l)^'Si 3 

72l72373lSl23'SfpS3P(73l)'^ y Sl2Sl3S3P(7l3)^ 2:3Sl3 

I23Z3S3P \ pNMHV 

• 513 ( 713 )^ 721'S13(73i)1 ) ^ 


(5.20) 


Sp[A^=4] (_p+; 1 -, 2 +, 3 -) = Sp(0) (-P+; ir, 2 + 3 -) 1/3 

732 2:2 S3P (723 ) ^ ( 731 ) ^ 7i3 T7NMHV 
2 ^2 
l23ZiSi3S23 

732 / 731^^2 _ 713721^^1^13 \ pNMHV 

(1 - 2;3 ) 723S3P V 721 73l5;3'Sl2(7l3)V ^ 

/ 712 ^ 13 ( 731 )^ ■§ 13 ( 731 )^ _ 2 : 1^13 _ 

V ^;3'S12S3P 5;3S123S3P 712732S123S3p(713)^ 


TTNMHV 

^3 


(5.21) 


{z3)[2P] ( 

[13] 


Sp[A^=4] (-P+;l+ 2-,3+) = Sp(0) (-P+;l+,2-,3+) (y, + (5.22) 

Sp[A^=4] (-P+;l-,2+,3+) = Sp(0) (-P+;lr,2+ 3+) {Vg + (5.23) 
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Sp[Rl (_p+; 1+, 2-, 3-) = Sp(°) (-P+; 1+ 2-, 3-) 



7 

2 



•512 


a2i{z2){z:i) ( _ 

[23] (12) y 272i73ia2isi23 


FbISc “ 7'Sl23'S23-f^2 (si2, S 123 ) 


1 


—3si23 — 


2^23 

732 




(5.24) 


Sp[Rl (-P+;ir,2+ 3 -) =Sp(°) (-P+;ir,2+,3-) (^) 

7 (5r\ ( 2 : 3 ) [ 2 P] ( /3l2/?2lSl3„im /3l2/32lsf3S23 


2 e V-■S12 


2 2 ] [13]si3 I 2731 S 123 


■nim 

■t^box ~ 


2731 S 123 


7^2 (Sl2, •S 123 ) 


+ 


{..., .1.3) + (»33, «3) 

2731 'SI 23 2732 

SI 2712 / 712 (2si 2S123 + 3S23S123 “ 3si2S23) 


732 


I ^ A ^23 

+ R 1 3712 - 

2 V 732^123 


2si2S123 

3 (712 - 1) S 12 
S 123 


+ ij 7^1 (523) S 123 ) 
— 1 ) Lq (•S 23 , ■S 123 ) 


7i 2'S23 , 712^23 , , (1 “ 712) 'SI2 (Sl23 “ S23) S12 

+ ;;-^ 712 + 


4732512 2732 S 123 


4si23 


4732 s' 


123 


(5.25) 


Sp[Rl (_p+; 1+,2-,3+) =Sp(o) (-P+;l+ 2-,3+) — 




€ V-512 


f I^R 
2 e 1—S12 


7 (5r a 


a 


3 

12 
3 

13 


^Im 

box 


“ 7 ( 65?2 + ^23 “ 6S12S123) L2 (Si2, S123) 


+ (3si23 — 2 S 23 ) Li (Si2, S 123 ) + -^Lq (si2, S 123 ) 

+ TkTfTnisli, (., 3 , ,,, 3 ) + (7^ - Til) i ,3 (, 33 ,,J 33 ) 


2c1i2 


a^2 CI 12 


(5123 + 523 )^ - 5i2 - 2S23512 _ 6^12 - ^23 ~ 2^23^12 513 _ 3 

4af2-®123523 4 si23512 0125123 2 


(5.26) 


Sp[Rl (-p+;ir,2+ 3+) =Sp(°) (-P+;ir,2+,3+) (^) 


l^-R 

2 e V-512 


7 (5r 012 / Fw CK?2S23 ? / N rv f / 

— R-H--1-— L 2 ( 512 , 5123 ) — 2012 S 23 F 1 (si 2 , 5123 ) 

2 2 013 y 2 2 

, 1 , (^123- 512 ) , 012523 of25i3 

+ 7-^0 (.512, 5i23j H- - -1-- 

2 2 si 23 5123 4 si 25 i 23 J J 


(5.27) 
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Spl'l (-P+;l±.2-,3-) = Sp<») (-P+; It,2,-. 3 -) 

- { - ^-^3 (si 2 , S123)- L2 (si 2 , S123) 

ci[12J(2ci) \ 723 721 


721731 ■S123 


-Ll (si 2 , S123) — 


721731'S123'S23 


-Lq (si2, S 123 ) + 


1 


2721'S12S123 


(5.28) 


Spm(-P+;l,-,2,+ ,3-)=Sp'»>(-P+;l,-,2+ 3 -) j 

+ ^ ^ 3 [ 13 ] -^ (2732S123L3 (si2, S123) — I2 (Sl2, S123) 


+ (^ + 2771? -A) 7 ( 3 .,, 3 ,,3) 

V732'S123 S12 512/ 

+ ( ^ U (.12, « 23 ) - 


I — ! — I 

V Si2'S123 'S12S123/ 


2si2S123 ) 


SpM (-P+;l+,2-,3+) =Sp(0) (-P+; 1+, 2“, 3+) 


10 

~9 


«12 

3®13 


(■^IS “ '^23) -^3 ('S12, 'S123)- L2 (Sl2, S123) 

ai2 


(5.29) 


3s 


H-(si2, S123) — Lq (si2, S123) — 7 


1 (•S123 — •513)^ — S13S123 — •S 12 (S 123 — S23) 


a 


12 


5125123 


513 ^13 (3si2 + 5si23 + 523) 


ai25l23 


2ai25l25l23 


(5.30) 


SpM (-P+; ir, 2+ 3+) = Sp(°) (-P+; 1^, 2+ 3+) 


10 

~9 


«12 / 
3ai3 y 

523 


— 2(112523^3 (512, 5i23) + 3ai2523l2 (512, 5123) + 2S13L1 (si2, S123) 


512 — 5123 


-Lq (si2, 5123 ) + 


5123 2si25i23^ 


523 «12523 


(5.31) 
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gp[scalar] ^ +^ 2“ , 3“) = Sp(°) (-P+; 1±, 2“ , 3“) 

Hz2){z^) 


^ 721732[12](23)si23^' (^12,^123) 


Spiscalar] (_p+; I" , 2+ 3“) = Sp® (-P+; Ir, 2+ 3“) ( 35 R- 2 ) 

(^ 3 ) ( , 3S2P im 3 si3S2P ; , . 3si3S2P ; , . 

+ :;-Fbox-Fi (si 2 , S 123 )-Fi (s23, ■S 123 ) 


[13](2P)l 2731S123 


731'S123 


731'S123 


_l_ 2si3g2P _ 2713732 Z2'S2P 27i2713^^2'S2P I 
732 S 12 S 123 723 (1 — 2 : 3 ) S 12 72323 S 12 j 


gp[scaiar] ^ 2“, 3+) = Sp(°) (-P+; 1^, 2", 3+) (^5n 


1 


at 


T-ilm 

^box 


(5.32) 


(5.33) 


'13 

+ SI 2 F 1 (s23, •S 123 ) + S 23 F 1 (si 2 , S 123 ) + Lo (si 2 , S 123 ) + Lq (s23, •5123)^ ^ (5.34) 


gp[scalar] ^_p+. ^-^ 2 + 3+) = Sp(°) (-P+; 1^, 2+ 3+) ^^(5r - 3ai2S23Ll (S 12 , S 123 ) ^ 

(5.35) 

The expressions for the non-zero independent one loop splitting amplitudes g —)■ qgq 
e 

Sp[A^=4] (-p+;l+ 2-,3-) =Sp(°) (-P+;1F,2-,3-) [ 


+ 


/ 721 ^ 23 ^ 23 ( 713 )^ 

V ZISIP 


1 ) - 


731'^123 (713«23)^ rpNMHV^ 

—Si—) 


(5.36) 


Sp[A^=4] (-P+;ir,2+,3+) = Sp(°) (-P+;ir,2+,3+) {Vg + 


(5.37) 
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Sp[scalar] ^_p+. 2 -^ 3 -) ^ Sp(0) ^ 3 -) 


2 

3~e 


10 

¥ 


+ ^ (si2, S123) - 1^0 (si 3 , S123) 

o 723 -J 



(5.38) 


— ■5123 


gp[scalar] ^_p+. ^ 2 +, 3+) = Sp(°) (-P+; , 2+, 3+) 

_R _|_ f 1^32 12 _ _!|-^ ^2 (Sl 2 , S 123 ) — 4ai2S23-£l ('Sl 2 , S 123 ) 

V ^ 7 i 2 ^ 7 i 3 / 


zy 

18 


+ 


s2 

*13 


S 123 


I If /„ „ N I 732*12 732 I 

+ 711^0 ('S12, Si 23 j + -r-2-"T^-r 7^- 

2 47f2Si23 712 47f3Si2Si23 2si2 


(5.39) 


The expressions for the others helicity configurations are obtained by conjugation operation 
or by reverting the fermion line, namely Sp (—T*; 1^, 2,3“^) = Sp y~P'^ 3^^, 2,1^^ . 


6 Super-symmetric Ward identities 

Super-symmetric Ward identities (SWI) allow us to relate amplitudes with different particle 
content in super-symmetric theories [67, 68]. We have checked that the super-symmetric 
primitive amplitudes presented in the previous sections obey the expected SWI which are 
easy to derive by taking the triple collinear limit of the equations presented [68]. The 
resulting equations differ for the MHV and NMHV helicity configurations. 


Sp[-^=4](_p+, 1+, 2-, 3+) - |^Sp[-^=^l(-P+, l4,2-, 3+) 

= 0, 

(6.1) 

Sp[-^=4](_p+, 1-, 2+, 3+) + |||sp[-^=^l(-P+, ir, 2+ 3+) 
Sp[-^=4](_p+, l+, 2-, 3-) - |^Sp[-^=^l(-P+, l4, 2-, 3-) 

= 0, 

(6.2) 

-|^Sp[-^=^l(-P+,l4,2-,3-) 

= 0, 

(6.3) 

Sp[-^=4](_p+, 1-, 2+, 3-) + |^Sp[-^=^l(-P+, ir, 2+ 3-) 



^ (|)spW=4](_p+,3-^2+,l-) 

= 0. 

(6.4) 


The SWI are also applicable to tree-level amplitudes in non-super-symmetric theories. 


7 Cross checks 

We check the universality of the splitting functions derived in the previous section numer¬ 
ically against the six parton amplitudes available in NJet [69]. In order to make sure 
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we could evaluate as close to the precise limit as possible, we implemented the checks in 
octuple precision using the QD and OneLoop [70] packages. 

We check the validity of Eq. (2.9) by computing the ratio between the two sides of the 
equation summed over the external helicities 


^collinearl23 — 


E l 

k=0 Z^k'=0 l^\p,X'p 






(7.1) 


where Ain,s and 'Pn,s in the denominator are defined from and Vn by summing over 
the external helicities: 




i(L,L') 

m-f 




{L,L') 
n—m+l,s\\p,\'p 


’(PM')) 


^ KHM-P’M))’ 

Al ,...,A 77 x 


Eq. (2.9) obviously implies 

^collinearl23 > t. I'.zj 

It is worth observing that the finite one-loop all-plus and all-minus four-gluon helicity 
amplitudes, while giving no contribution to the NLO squared matrix element, they give 
instead a finite contribution to rcoiiineari 23 because of spin correlations. 

In Fig. 5 we plot rcoiiineari 23 — 1 as a function of the invariant mass S 123 of the three 
collinear partons. More in detail we verify the validity of Eq. (7.2) in double, double¬ 
double and double-quadruple precision for both gluon (on the left) and quark (on the right) 
splitting functions. As one can see, going to higher precision allowed us to make stronger 
checks on phase-space space points which are closer to the limit, where the numerical 
evaluation is highly unstable at lower floating-point precision. 

Similarly, we also numerically verified Eq. (2.8) for each primitive amplitude and all 
the helicity conhgurations, although all of these already contribute to the check described 
above. 

As well as the numerical checks we have also verified that all splitting functions factorise 
correctly in the iterated collinear limit, 


limSp^-^), 1^1,2^2^ 3A3) 

L 

= EE Sp(^-^), 3^3)Sp('^)(-Pf/, 2 ^ 2 ), (7.3) 

A=± fc=0 


where the scale S 12 S 123 and = P 12 +P 3 . All di-logarithms drop out in this limit 

though some care should be taken to ensure the hierarchy of scales is imposed correctly. 
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Figure 5: Numerical check of the collinear limit of rcoiiineari 23 — 1, with rcoiiineari 23 defined 
by Eq. (7.1), as a function of the invariant mass S 123 , from S 123 = 10“^ to S 123 = 10“^°. 
For this check we set n/ = 5, y/s = 10^ and = 10^/7. The plot on the left shows the 
all-gluon case, while the one on the right shows the quark case. In the latter it was not 
possible to obtain numerical results in double precision for S 123 < 10 “^®. 


8 Conclusions 

In this article we have computed the one-loop triple collinear splitting amplitudes in QCD 
initiated by a gluon. These functions are one of the last remaining ingredients to complete 
the classification of universal infrared limits relevant at N^LO. 

Some effort has been taken to ensure the splitting amplitudes have compact analytic 
forms. We made use of the spinor-helicity formalism and super-symmetric decompositions 
and related the pure gluonic amplitudes to the ones containing a quark anti-quark pair. The 
primitive amplitude colour decomposition was also a useful tool to express full colour and 
helicity summed splitting functions which were all checked explicitly against the numerical 
matrix elements for 2 —)> 4 scattering in NJet. In the course of these checks we made 
use of the high precision numerical evaluation available with up to 64 digits via the QD 
package. This allowed us to probe deep into the collinear limit and verify that all parts of 
the computation behaved correctly. This was particularly important for the spin correlated 
and sub-leading colour corrections which are significantly suppressed. 

There are still some missing ingredients needed for the constructions of a fully differ¬ 
ential N^LO subtraction scheme. Firstly, the quark initiated channels are still unavailable - 
they are not directly accessible from the H + 2j amplitudes since they have been computed 
in the effective theory where the Higgs couples only to gluons. The necessary splitting 
amplitudes could be extracted from the vector boson plus four parton one-loop amplitudes 
[71]. 

Secondly when integrating the splitting functions over the unresolved phase space the 
expansion of the limit may be required to higher order in the dimensional regularisation 
parameter e. This would require a new computation of the one-loop matrix elements valid 
in 77 = 4 — 2e dimensions which is quite feasible using modern unitarity methods. The 
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appearance of the one-loop pentagon function in the full D-dimensional amplitude may 
complicate this part of the computation even if it is only required in the triple collinear 
limit. 

We hope that the expressions presented here will be of use in future high precision 
QCD computations. 
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A Generation of collinear phase space points 


In this Appedix we illustrate a practical way to generate a set of on-shell re-particle phase- 
space points where the first m particles approach the collinear limit 1|| • • • ||rre. The limit is 
approached by varying a single free parameter d as 5 —)• 0 and it is based on the parametri- 
sation presented in Section 3. This has been used for the numerical checks we discussed in 
Section 7. 

As a first step we generate an on-shell (re — m -|- l)-particle phase space point defining 
the set of momenta 

{P, Pm+l(0), Pm-H2(0), . . . , Pn(0)} (A.l) 

where, as suggested by the notation, pi{0) for i > m + 1 are the momenta of the non- 
collinear particles at d = 0, while P is the sum of the collinear momenta in the limit. We 
then define the exact collinear limit as the set of momenta 


{ziP, Z 2 p, . . . , ZmP, Pm+pO), • ■ ■ , Pn(0)}, (A.2) 


where zt are randomly generated real numbers satisfying Eq. (3.6). In order to avoid regions 
with soft kinematics (which would introduce other kinds of singularities) one can generate 
a set of random numbers between, for example, 1 and 3 and divide them by their sum. 

In order to define the orthogonal direction we must specify the reference vector r] 
appearing in Eq. (3.2). A particular convenient choice is one of the non-collinear vectors, 
i.e. 

r]^=P^+i{0). (A.3) 

The orthogonal direction is thus spanned by the two complex vectors 

{Ph^lv] {vh^\P] ., ,X 

2 ’ 2 ■ ^ ’ 

While these are particularly convenient when working with the spinor-helicity formalism, 
for numerical checks with real kinematics it is convenient to define two real linear combi¬ 
nations 


^ _ 1 ( jPppr^] , {rj\jpP] \ ^ _ 1 ( {Phpy] (r/|7^|P] ^ 

2 y 2 2 j ' 2 2 j ' 


(A.5) 
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Hence the orthogonal vectors kit - are defined as 

1 ,Z 


= Vki + 1/2,i vl^,± (A.6) 

where yi^i and 2 / 2 ,i are randomly generated real numbers satisfying 

X] 2/bi = X] ^2,* = 0. (A.7) 

i i 

The variables 2 / 1 ,i and 2 / 2,1 are related to the spinor variables (zj), [zj], (cj,) and [ujj\ intro¬ 
duced in Section 3 by 

2/i,i = {zi) [uJi] + {uji) [zi], 2 / 2 ,i = i {{zi) [uJi] - (uji) [zi]), (A.8) 

as one can check by requiring consistency with Eq. (3.9). As already stated, these spinor 
variables differ by a phase from the usual parametrisation in terms of If 

{zi) = [ziY = {uJi) = [wi]* = y/uTie^'l’ (A.9) 

then 

yi^i = 2y/zYuHcos{<j)-6), 2 / 2 ,i = 2 y/zi oji sm{(j) - 0). (A.10) 

From here it is easy to see that while the parametrisation in terms of {zi) and [z*] has the 
advantage of producing results that are analytic functions of the spinor variables in the 
complex plane, the parametrisation in terms of is in fact entirely equivalent in the 
physical region. 

Using j as in Eq. (A.6), one can simply define the momenta pi, ■ ■ ■ ,Pm for any value 

of the free parameter 5 using Eq. (3.2). With our choice of p we can absorb the recoil by 

defining 

m ^2 ^2 

P^+l{S) = (1 + Z-7^^) /5m+l(0) (A.ll) 

p^{5)=p^{0), m + 2 <i <n. (A.12) 
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